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ABSTRACT 



The established concept of temperature determination from dual spectral band 
radiometric measurement is analyzed for the case that the measurement device is sensi- 
tive to total photon emittance rather than radiant cmittance. A temperature dependent 
ratio of measurements is calculated for two distinct spectral bands using a black body 
distribution. It is shown that, if the bands are not self contained, then there exists a 
one-to-one correspondence between temperature and ratio. 

A prototype algorithm is proposed and tested which demonstrates the feasibility of 
extracting both the temperature and range by using three distinct spectral band meas- 
urements. The model assumes nonhomogeneous, wavelength-dependent atmospheric 
attenuation. The target and sensor systems each have an arbitrary location in a vertical 
plane relative to the earth. In the computer simulations, LOWTRAN data was used. 
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I. INTRODUCTION 



A. REVIEW OF PRIOR WORK 

In the photon detection process electromagnetic radiation causes three kinds of 
mechanism such as photon effects, thermal effects, and wave interaction effects 
[Ref. 1]. The focus of the remaining discussion will be confined to the first two effects 
which cover the majority of applications involving the detection of thermal radiation. 
One common photon effect is photoconductivity. Photoconductivity is observed in 
semiconductors and characterized by photon energy and energy gap. One common de- 
vice which uses a thermal effect is the bolometer. Bolometers can be made of any ma- 
terial which has a temperature-dependent resistance. This effect is characterized by the 
change of resistance from the heating effect of incident radiation. Therefore, Plank's 
equation for radiant emittance can only be applied directly to thermal effects. Although 
some devices have a detection mechanism based on the thermal heating effects due to 
radiation, this is not true of semiconductor-based detectors. Semiconductors devices are 
characterized by an energy gap. Only photons with energies greater than this gap can 
produce electron hole pairs which contribute to the detected signal. The excess photon 
energy is dissipated as heat [Ref 2] and therefore does not contribute to the electrical 
signal generated. Specifically the integrated spectral emittance [Ref 3], i.e., photon flux, 
is a more appropriate signal indicator than the radiant emittance for measurements 
based on semiconductor devices. 

R.B. Johnson and E.E. Branstetter presented a numerical method for the integration 
of Plank's equation [Ref 4]. Johnson derived an approximated equation for effective 
radiant sterance [Ref 5]. Normalization of spectral bandwidth, for convenience in 
analysis, was provided by F.E. Nicodemus [Ref 6]. Based on Plank's equation for ra- 
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diant emittance, M.II. Horman presented the technique for temperature determination 
using the ratio of output from two bands (Ref. 7]. On the other hand, analysis of the 
dual spectral band temperature measurement for devices based on the photon effects has 
not been developed. 

One standard method for the determination of the gray-body emissivity and the 
temperature of an object is based on two separate narrowband radiometric power 
measurements. In combination with Plank's radiation equation for a gray-body, two 
narrowband measurements are sufiicient to produce a unique value for both temperature 
and emissivity. For practical reasons, such as improvement of signal-to-noise ratio, 
most instrumentation take measurements over bands which are too wide to be consid- 
ered narrowband. In this more general situation the question of uniqueness of temper- 
ature requires more careful consideration. This issue has been recently analyzed for the 
distribution based on a Plank's spectral radiant emittance [Ref 8]. The authors devel- 
oped an analysis based on the evaluation of the ratio (^) of two measurements of the 
integrated power distribution over separate bands. Three cases, nonoverlapping bands, 
partially overlapping bands and totally overlapping bands, were considered. The analyt- 
ical investigation demonstrated that the last case does not always generate a one-to-one 
correspondence between the ratio ^ and the temperature. 

For the passive ranging technique, P.J. Ovrevo and R.C. Wood presented the 
method using absorption property of the atmosphere [Ref 9], and J.R. Jenness, Jr., and 
F.J. Shimukonis proposed the principle of passive ranging using the ratio of the signals 
in narrow spectral bands received [Ref 10]. 

For optical properties of the atmosphere, R.A. McClatchey and others described 
effects of index of refraction in the atmosphere [Ref 11: pp. 4Ij. It was found that 
transmission factors are affected by no more than 1 - 2 % if refraction is included in the 
calculation. 
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B. OVERVIEW OF THE THESIS 



Chapter I is the introduction and covers both a review of relevant literature and a 
thesis overview. In the interest of keeping this thesis fairly self-contained, Chapter II 
provides a brief technical description of essential background preUminaries. In Chapter 
III, effects due to atmospheric attenuation and detector responsivity are ignored. Using 
dual spectral measurements, it is shown that it is possible to uniquely determine a tar- 
get's temperature if the detector bands are not 100 % overlapping. Specifically, the ratio 
of the dual spectral band measurements is shown to exhibit a mathematically predictable 
one-to-one correspondence with the temperature, as long as the bands are not totally 
contained within the other. 

In Chapter IV, a numerical program for calculating a spectral band measurement 
without the assumption of constant atmospheric attenuation and detector responsivity 
is discussed and tested In the process of this development, a novel definition for a 
"responsivity", applicable to devices sensitive to photon emittance, has been introduced. 
If the effects of atmospheric attenuation on the dual spectral band measurement are not 
equivalent, then the ratio scheme for temperature determination, discussed in Chapter 
III, is not applicable. It is precisely under these conditions that by using three distinct 
bands both the temperature and the range can be determined. 

In Chapter V, a prototype algorithm is discussed, which demonstrates the feasibility 
of obtaining both the temperature and the range of a target using three distinct spectral 
band measurements. Complications introduced by not assuming constant altitude of the 
target-to-sensor trajectory are addressed. LOWTRAN data was also used here to char- 
acterize atmospheric attenuation. Conclusions for the thesis are presented in Chapter 
VI. In order to not obscure the main points of the thesis, the extensive mathematical 
derivations have been relegated to appendices. Both Appendix A and Appendix B fit this 
description. Appendix C, Appendix D and Appendix E have been created for those 
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readers who are interested in an algorithmic description of programs used in the com- 
puter simulation. 
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11. BACKGROUND 



In this chapter, an alternate approach based on photon emittance rather than radi- 
ant eniittance is applied. This kind of treatment is appropriate for photovaltaic or 
photoconductive semiconductor devices but not radiometers. First, a few basic as- 
sumptions need to be presented. The discussion starts from the ideal case with the fol- 
lowing preliminary assumptions [Ref 8: pp. 1256]: 

• opaque object 

• Lambertian surface 

• object in thermal equilibrium 

• spectrally constant emissivity (gray body) 

• photon emittance derived from Plank's radiation law 

• unity spectral transmittance of intervening media between the object and measuring 
instrument. 

Plank's radiation law gives the spectral radiant emittance, which is given by 
[Ref 3; pp. 35]: 



W,{X,T) 



2nt{X)hc^ 

X\exp{hclXkl')-\] 



^2 ] 

[IV • cm • jxm ], 



( 2 . 1 ) 



where c(A) is the spectral emissivity, h is Plank's constant, X is the wavelength, c is the 
speed of light, k is Boltzmann's constant, and T is the object's temperature. In this sec- 
tion, since the emissivity is assumed to be spectrally constant, c(2) is denoted as e. 

The photon emittance is derived from (2.1) dividing by hejX, which is the energy 
associated with one photon [Ref 3: pp. 38]. Therefore photon emittance is the flux re- 
lated to the number of photons. The resultant photon emittance equation is obtained 
as 
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2ncc 



\photons • s ' • cm 



X\exp{hcjXkT) - 1] 



-2 



•” 1 1 

pm ] 



( 2 . 2 ) 



The integral of (2.2) over some band limits gives the photon emittance in the bands. 
The photon emittance within the spectral band limits [/,, X^ is given by 



QiX2,X„T) = j\{X,T)dX (2.3) 

From the assumptions, emissivity is assumed to be spectrally constant. Therefore 
photon emittance for a gray body is simply expressed as photon emittance for a black 
body multiplied by a constant emissivity and given by 

Qa ~ ^QbBx' ( 2 ' 4 ) 



From (2.3), it follows from (2.4) that 



Q{X2, X],T) — 1[ 

^2 (2.5) 

= ^Qbb{^2^ 

where >(2 < and T) is the spectral photon emittance for a blackbody [e(.^) =1]. 

It has been shown that the temperature of an object can be determined by taking a 
ratio of measured radiant emittance of an object in each of two spectral bands. 
[Ref 8: pp. 1256] In Chapter II and Chapter III, this concept is applied in a similar 
manner except the assumption is made that a ratio is formed from the measured photon 
emittance. Let the two spectral bands be denoted as [X 2 , ij and [A4, /Ij] . Let the ratio 
for photon emittance be defined as 
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( 2 . 6 ) 



■^ 3 » ■^ 4 > ~ 



QkKhJ) ■ 



Since from the assumption, emissivity is uniform over the spectral band, emissivity is 
cancelled out, and the ratio becomes the ratio of the photon emittances associated with 
a black body, as shovm by 



■^ 2 » ^ 3 » ■^•41 ^ 



QbbK^-1^1 



(2.7) 



This implies that is a function of the object's temperature T only and the ratio can 
indicate the object temperature. When the dual spectral-band photon emittance method 
is used, it is critical to know that the ratio and the temperature have a one-to-one 
correspondence. For this purposes, three general cases of spectral bands which cover 

all possibilities are considered. The possible cases are [Ref 8: pp. 1256] 

1. non-overlapping bands (Figure la), 

2. partially overlapping bands and common band limits (Figure lb,c), and 

3. totally contained band (Figure Id). 

These three cases are presented graphically in Figure 1. In the following two chapters 
the validity of temperature determination based on dual spectral photon emittance for 
the ideal case is discussed. 
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a. nonoverlapping case 




b. partially overlapping case 



13 2 4 

c. one common band limit case 



1,3 2 4 

d. totally contained case 



13 4 2 



Figure 1. graphical representation of possible band combination: three cases are 
nonoverlapping, partially overlapping, and totally contained. 
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III. DUAL SPECTRAL BAND MEASUREMENTS FOR TEMPERATURE 



DETERMINATION 

In this section the ratio of photon emittance defined in (2.6) is considered for the 
various cases discussed at the end of Chapter II. In order to validate the dual spectral 
bands method, it is required to show one-to-one correspondence between the ratio and 
the temperature of the object. For this purpose, the derivative of the ratio of band lim- 
ited photon fluxes, in (2.7), with respect to temperature is performed. If the ratio of 

photon emittance is proved to be either strictly increasing or strictly decreasing 

function with respect to temperature, i.e., either dSiJdT > 0 or d^JdT < 0 for all T, then 
the ratio of photon emittances and the temperature exhibit a one-to-one correspond- 
ence. It follows that this approach is potentially applicable to the temperature determi- 
nation of grey bodies. 

To simplify the integral of (2.2) , the variables of integration in photon emittance 
Aj,T) are replaced by the dimensionless parameters [Ref 8: pp. 1257] 

and 

a=-^, (3.2a) 

where T, is an arbitrary reference temperature. By prior convention [Ref 8: pp. 1257], 
it is chosen to be 



Tr = 



he 

A^k ’ 



{3.2b) 
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where yl„ as introduced in Chapter II, is the highest wavelength in spectral band 1. Using 
(3.2b) as a reference temperature in (3.1), the spectral band limits of ij/ are given by: 



he 



‘ ^-ikTr 



X, = 



(3.3) 



where .1, is one of the spectral band limits [Ref. 8: pp. 1257]. After substitutions into 
(2.5), the new representation for the photon emittance is given by 






(3.4) 



where [Ref 8: pp. 1257] 



gW = 



r 



exp(iA) - 1 



(3.5) 



The ratio of photon emittance is obtained as 



0tp = 



f«2 

J. * 






faJC4 

J g{\l/)d\j/ 






(3.6) 



and the derivative of the ratio of photon emittance signals (3.6) with respect to temper- 
ature comes out to be 



J 'oUj foLXt CaX2 faj(4 

gWd^ - gWd^ A'l')d^ 



dT T 


r 


o 




gWdiP 


z 




L 
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where 



cxp{\j/)\l> 
[exp{iP)-lf ■ 



(The derivation of (3.7) is in Appendix A.) To simplify the expression of (3.7), the fol- 
lowing operator is defined [Ref 8: pp. 1257]. 



[a,^)](l, X 2 , ^ 3 , ^ 4 , a) = f a(iA)# f °^6(^)# - f '^^(^)# a{\{/)d\l/ (3.9) 

*^a *^a *'ajcj 



Then the ratio becomes 






BT T 


f<W4 








2 




L‘'«3 J 





(3.10) 



Since the denominator of (3.10) is always positive, the sign of d^JdT is determined by 
the sign of the operator [f,gj. It is necessary to show that the sign of [f,g], which is the 
sign of (3.10), does not change with a, or equivalently with temperature. The temper- 
ature dependence of [f,g] appears only in a, which can be seen in the limits of integration. 
The integrands do not have an explicit dependence on temperature. 

A. NON-OVERLAPPING CASE 

In order to prove that the sign of [f,g] is fixed over the domain of a, several proper- 
ties o[J[\l/) andg(^) are required. They are [Ref 8; pp. 1257] 

1. yio)=^(o)=o 

2. 0 < <y(iA) for 4/>0 

3. = linV-oog(<A) = 0 

4. For p> l,y(iA) — figi'P) has exactly one root on the interval 4f e (0, oo). 
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The discussion to follow will clarify the significance of these requirements. Grapliical 
representation is useful to clarify these properties and it is shown in Figure 2. This figure 
is especially useful for visualization of the fourth property. For values of > 1, it is ob- 
vious that crosses Pg{i}/) exactly once. A more formal discussion of the four prop- 
erties is provided in the next paragraph. 

The first property is proven by using THopital's rule. In order to prove the second 
property, define 






Ml 



(3.11a) 



which by direct substitution from (3.5) and (3.8) is given by 



PW 



1 - exp( - \]/) ■ 



(3. nil) 



It then can be shown that P{\1/) > 1 for > 0. This follows from the derivative of P(^), 
which is : 



dP{ip) 1 — exp( — ij/) — exp( — \p) 
# ~ (l-exp(-^^))^ 



After expanding (3.12a), it is clear that ; 

dPW _ 1 ^ exp( -\J/) ^ ^ 

dll/ l-expi-ip) (1 _ exp( - i/i))^ 



(3.125) 



which follows from the observation that the second term in the expansion is smaller in 
magnitude. Also, by direct application of rHopital's rule on (3.11b), lim^^o+P(i/i) = 1. 
These two results imply that F(^) > 1. Thusy((/^) > giip). Consequently, after noting that 
bothy(^), as given in (3.8), and g{ip), as given in (3.5), are strictly positive, property two 
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follows. The third property is simply proven by limit evaluation. For the fourth prop- 
erty, it is necessary to evaluate J[\p) — Pg{\]/) given by 



A^) - PgW = - — 77; — “T - P — T 

( exp lA) - 1) - * 



(3.13a) 



which can be expressed as 



A'l') - PgW 






/ * 1 ' 

exp(iA) - 1 V exp(iA) - 1 



p 



(3.13A) 



By noting (3.1 lb), it follows that 

I (3.13c) 

Since lin\^o*!2('A) = 1 and p is constrained to be greater than one, the quantity 

A'P) ~ Pgi'P) satisfies, after reference to (3.1 la), the inequality 

A'l')-PgW<0 (small ^) (3.14a) 

for sufficiently small \p. Also since lim^^oooF(i/() = 00 , the inequality has to reverse for 
large enough \p. In other words 

A^)-PgW>0 (large lA). (3.14Zi) 

Since, from (3.12a), P{\p) is a strictly increasing function, it is correct to conclude the 
inequality reversal can only occur once. Property 4 is therefore proven. 

For the non- overlapping case represented on Figure la, there are only two sets of 
relations which can define the spectral bands (3.3): 

\<X 2 ^X 2 <x^, (3.15a) 
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and 



^ 1 < Jfj* (3. 156) 

Relation (3.15a) is considered first. Motivated by the fourth property previously stated, 
let 



/?(«) = 



A°^X2) 
g{<^2) ’ 



(3.16) 



The quantities y(i/^) and p{a.)g{\l/) must cross each other only once at the value of \jj de- 
noted as X 2 - This is represented in Figure 3. The geometrical aspects of figures dictates 
the following relations for all values of a; 







(3.17a) 



and 




00C3 



(3.176) 



Substitution (3.17a) and (3.17b) into (3.9) shows the operator [f,g] satisfies: 



CaXj Cixxa f“^ 

X 2 , X 3 , a) < p{a)g{\l/)chj/ g{\p)di^ - g(i/')^^ p{a)g{\l/)d\p = 0 

‘'oJCj •'a *^aX3 



(3.18) 
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for all the values of a (3.2a). According to the brief description following (3.10), this 
proves that d^JdT <0 for all the values of T and therefore is a strictly decreasing 
function of temperature. For the second relation of non- overlapping bands (3.15b), a 
similar analysis would show that d^JdT > 0 and is a strictly increasing function of 
temperature. 

B. PARTIALLY OVERLAPPING AND COMMON BAND LIMIT CASE 

In this section, the case which is partially overlapping and the common band limit 
case are considered. Two sets of relations are considered in this case and they are 

1 Jf3 < Jf2 < ^4 or 1 < < Jf2 ^ -^4 (3. 19<j) 

X3 < 1 < ^4 < ac2 or jr3 < 1 < jt 4 ^ ^2 (3.19f>) 

For this discussion, the graphical method previously introduced is useful [Ref. 8: pp. 
1258]. Figure 4 illustrates the first set of the relations (3.19a) for partial overlap. With 
reference to Figure 4, six sets of integrals are defined as follows: 







f<«3 



{3.20a, b) 



S2 = 




f<«2 

^2 = J f('A)# 

axj 



{3.20c, d) 



rax* raX4 


{3.20eJ) 


^olX2 





After inserting this formalism into (3.9), the bracketed operator [f,g] becomes 



15 



[f,g]il,X 2 , X 3 , x^, a) = (S, + + ^ 3 ) ~ + ^2)(*52 + ^a)- 



(3.21a) 



After expansion and recombination of (3.21a), it follows 



Ml. ^ 2 . ^ 3 . ^4. a) = + S 2 ) - + ^ 2 )] + (Si^2 - ^i^i). (3.216) 

This can be retranslated using (3.20) into two terms which follow 

l/lg](l,X 2 , Xj, X 4 , a) = [/;g](l, X 2 , X 2 , X 4 , a) + l/lg](l, X 3 , xj, X 2 , a). (3.21c) 

Case (3.19a) will now be applied on expression (3.21c). If 1 < jcj < jfj < jr^, then both 
terms in (3.21c) are of the form 1 <X 2 = Xj<x^ which satisfies the nonoverlap condition 
(3.15a). Since, as indicated in the previous section, the particular nonoverlap condition 
(3.15a) generates negative values for [f,g], the expression (3.21c) must be negative for the 
stated inequality. 

For common band limit case, which is 1 = jtj < jtj < or 1 < 0:3 < Xj = 3 ^ 4 , one of the 
terms in ( 3 . 21 c) becomes zero. However, the other term remains negative, and therefore 
[f,g] remains negative. This discussion implies that the ratio is a strictly decreasing 
function of temperature. Similarly, for the second relation of this case (3.19b), the ratio 
is a strictly increasing function of temperature. 

C. TOTALLY CONTAINED SPECTRAL-BAND CASE 

This section discusses case 3, where one interval is entirely contained within the 
other. By the use of asymptotic approximations (Ref. 8 ; pp. 1258), it is possible to show 
that the ratio may have the same value for two or more temperatures. Therefore, the 
measurement to determine the temperature is not unique. Two examples, which are the 
case where a is large and the case where a is small, are discussed. 

As previously discussed, the sign of d^JdT is determined by the sign of [f,g]. There 
are two sets of relations which are 1 < oTj < jt 4 < jTj ai^d jCj < 1 < jTj < x^. The first set of 
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relations is arbitrarily chosen for discussion since the approach is similar for the other 
relation. For large or small value of a asymptotic methods are useful to determine the 
sign of [f,g]. (The derivation using asymptotic approximation is in Appendix B.) For the 
large value of a (i.e., small temperature ), [f,g] is approximated as 

[/;g]2::a^jf3(l - .X3) expC - a(jf3 + 1)]. (3.22) 

Since jTj is greater than one, [f,g] becomes negative. For small values of a (i.e., large 
temperature ), [f,g] is approximated as: 

^ - x^3){xl - 1) - {xl -\){xl - ^3')]. (3.23) 

Given and x^, it is always possible to find a sufliciently large such that (3.23) is 
positive. Combining the tw'o results for small a and large a, the sign of dSiJdT is not 
fixed over a (i.e., temperature ). This means that the ratio does not have a one-to-one 
correspondence with temperature. Therefore, for the case when one spectral band is 
completely contained within the other, the temperature determination method based on 
the dual spectral band measurement of photon emittance is not, in general, valid. 

D. RATIO CALCULATION FOR PHOTON EMITTANCE 

In order to generate the output ratio calculation, the numerical method for inte- 
gration known as Simpson's rule [Ref 12: pp. 95] was used. The flow chart for this 
calculation is shown in Appendix C. The examples were demonstrated using the condi- 
tions in [Ref 8: pp. 1258-1259). One case has a band combination of [2.5 5.5 /mi] 

and [3.9 /mi, 4.1 /im| shown in Figure 5. The second case has a band combination of 
[10.0 /an, 14.4 /im] and [11.8 /im, 12.2 /im] shown in Figure 6. In these figures, the nor- 
malized ratio denoted as was used to enhance the results. Normalization is done by 
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(3.24) 



~ ^min 

•^P ~ ~3 , 

“^max *^mm 

where is the maximum value of and is the minimum value of Both cases 
are slightly different from the results in [Ref 8: pp 1259], but the general features are the 
same. One ratio value may correspond to two or more temperature values as predicted 
in Section C. These figures indicate that the minima of the ratio occur in different po- 
sitions depending on the bands selected. 

E. DISCUSSION OF COMPUTER RESULTS 

Computer simulations were performed for the three cases discussed in this chapter. 
The first example is the nonoverlapping case and the band pair is chosen to be [3.0 ^m, 
4.0 /xm] and [4.0 /xm, 5.0 /xm]. Figure 7 shows this example and it exhibits the decreasing 
feature as expected in Section A. The second example is the partially overlapping case 
and the band pair is chosen to be [3.0 /xm, 4.5 /xm] and [3.5 /xm, 5.0 /xm]. Figure 8 shows 
this case and the results agree with the discussion in Section B. The third example is the 
case of totally contained band within the other and the band pair is [3.0 /xm, 5.0 /xm] and 
[3.5 /xm, 4.5 /xm]. Figure 9 represents this case and the results show that one-to-one 
correspondence does not exist in this case. It can be seen from Figure 9 that a ratio value 
of 2.0 corresponds to both 460 degrees Kelvin and 1400 degrees. 
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Figure 2. Plot of f(^) and Pg(^) for U2y...,5 
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Figure 3. Nonoverlapping case.: The constant ^(a) is chosen such that 
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Figure 4. Partially overlapping case or common band limit case: Each region is 

defined in the text. 
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Figure 5. Example of totally contained case (normalized ratio): Bands are [2.5 

^iin, 5.5 and [3.9 ^lm, 4.1 ^m]. 
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Figure 6. Example of totally contained case (normalized ratio): Bands are [10 
/im, 14.4 /im] and [11.8 tim, 12.2 ^m]. 
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Figure 7. Example of the non-overlapping case: Bands are [4.0 nm, 5.0 Mm] and 
[3.0 Mm, 4.0 Mm]. 
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Figure 8. Example of the partially overlapping case: Bands are [3.5 /zm, 5.0 tim] 

and [3.0 4.5 /xnr]. 
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Figure 9. Example of the totally contained case: Bands are [3.0 ^m, 5.0 ^m] and 

[3.5 ^m, 4.5 /ini]. 
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IV. A GENERAL MODEL FOR TEMPERATURE DETERMINATION VIA 



DUAL SPECTRAL BAND MEASUREMENTS 

In Chapter II, ideal photon emittance temperature determination was discussed. 
Since this discussion was dealing with the ratio without considering measurement of the 
photon emittance as a system, it was not always practical. Therefore, some other factors 
should be taken into account to realize this method. They are emissivity of the object, 
atmospheric transmittance, detector responsivity, and so on. In the previous discussion, 
those parameters were assumed to be spectrally constant or unity. Strictly speaking, 
they depend on the wavelength. Also after considering these factors as part of the sys- 
tem, the resultant ratio becomes the ratio of the output voltages from the detectors of 
two bands instead of the ratio of photon emittance. A brief discussion of these param- 
eters follow. 

A. EMISSIVITY 

The emissivity dictates how much thermal energy is radiated from an object. Its 
value lies between zero and unity depending upon the material characteristics of thermal 
source. It is denoted as e(/ 1). Although emissivity for some materials has temperature 
dependency, it is assumed to be independent of temperature to make discussion simpli- 
fied. Materials are divided into three types based on the value of emissivity. They are 

[Ref. 3: pp. 40] 

1. A black body, e(yi) = 1. 

2. A gray body, c(2) = constant < 1. 

3. A selective absorber, e(2) varies with wavelength. 

The spectrally constant emissivity cases were covered in the previous Chapters II and 
III. Therefore, wavelength-dependent emissivity is the primary focus of discussion in this 
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chapter. In addition, since emissivity depends on the types of materials, it is not possible 
to get reliable analytical expressions for all selective absorbers. Usually emissivity values 
are obtained by measurement. 

B. ATMOSPHERIC TRANSMITTANCE 

The atmosphere is known to behave as a wavelength-selective absorber. This im- 
plies that the transmission of signals through the atmosphere becomes wavelength de- 
pendent. Specifically, the atmosphere affects the transmission of infrared radiation. The 
main compositions of the atmosphere are nitrogen, carbon dioxide, oxygen, and water 
vapor. Among these, molecules which have nonsymmetric molecular structure are pri- 
mary responsible for absorption. Thus carbon dioxide and water vapor are the dominant 
sources for absorption. In addition some other small constituents of the atmosphere, 
which have non-symmetric structure, also contribute absorption. 

Atmospheric scattering is the other important phenomenon which affects atmo- 
spheric transmission. The strongest atmospheric scattering occurs in the case when the 
size of aerosols and the wavelength of the signal are almost the same. In this particular 
case the scattering coefficient depends highly on a w’avelength of the signal. The Mie 
theory' describes this phenomenon. 

The atmospheric transmittance is determined by atmospheric absorption coefficient 
and atmospheric scattering coefficient. It is given by 

T(^;c) = exp[-(a(.l)-fy(^)W (4.1) 

where a(>l) is atmospheric absorption coefficient, y(yl) is atmospheric scattering coeffi- 
cient, and X is the distance betw’een the object and the detector. An example 
transmittance data by the NFS LOWTRAN6 program [Ref 13] at 20 meters in height, 
0 degree elevation for 1 km and 5 km distances under 1962 U.S. standard atmospheric 
conditions are shown in Figure 10. The importance here is that the transmittance de- 
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pends not only the wavelength but also distance and atmospheric conditions. Moreover, 
atmospheric conditions vary' with height and thus the transmission depends on the path 
between the object and the detector. Unless the path and atmospheric conditions are 
known, it is impossible to apply the ideal ratio calculation. Therefore, the transmission 
path and atmospheric conditions have to be provided in order to obtain the atmospheric 
transmittance. A recent report [Ref. 11: pp. 41] on transmittance calculation for the at- 
mosphere has concluded that refraction effects can essentially be ignored. 

C. DETECTOR RESPONSIVITY 

The detector responsivity is generally obtained by dividing the measured output 
signal voltage by the optical power falling on the detector. And it is defined as 



R = 



Kn,s 

HA, 



(4.2) 



where R is the responsivity [V/W], is the RMS output voltage, H is the RMS 
irradiance [Wlcm^] at the detector location, and A, is the detector area [cm^]. However, 
this quantity is not the appropriate parameter for the development of the photon 
emittance case. A new proposed responsivity R' is defined by modifying (4.2) in order 
to adjust it for the photon emittance calculation. It is given by dividing the measured 
output voltage signal by the number of photons into the detector. The proposed R' is 
defined as 



R' = 






(4.3) 



where R.' is the photon responsivity [ V • photons^' • 5 ], and Q is the photon emittance 
[photons • s~^ • cnr'^\. It follow's from definitions (4.2) and (4.3) that the photon 
responsivity R' is simply obtained by multiplying responsivity R by hcjA, 
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/?' = /? -y. 



(4.4) 



as shown above. For photon detectors sensitive to the number of photons, R' should 
be fairly constant out to the cut off wavelength. 

Unfortunately, detector responsivity (4.2) is not usually provided among detector 
specifications, but the specific detectivity D* is given instead. The D* is the normalized 
quantity of detectivity defined by [Ref. 3: pp. 270]: 



where B is the electrical system bandwidth in heltz, is the area of the detector, and 
D is the detector detectivity. Tliis quantity is often expressed in terms of the reciprocal 
of the noise equivalent power (NEP). Therefore, the D* is also given by [Ref 3: pp. 
270] 



The noise equivalent power is the radiant flux required to obtain an output signal equal 
to the noise of the detector. It is given by 



where V, is RMS signal voltage and is RMS noise voltage. From (4.6), the specific 
detectivity becomes 



D* = DJA^D 



(4.5u) 




(4.6) 




(4.7) 
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The photon responsivity is obtained by combining (4.3), (4.4), and (4.7), and it is 



D*Vn he 

JaJ -1 ' 



(4.8) 



This equation shows the specific detectivity is closely related to the photon responsivity. 
From the expected wavelength dependencies in (4.8), it follows that the quantity D* 
should increase linearly with wavelength out to the cutoff wavelength [Ref 3: pp. 294]. 
Since photon detectors respond to the number of photons absorbed in the detector up 
to cutoff wavelength (usually dictated by the band-gap), the modified responsivity of 
photon detectors is ideally independent of wavelength. 

D. COMPUTER SIMULATION OF TEMPERATURE DETERMINATION 

Based on above discussion, the ratio of photon emittance for black body is replaced 
by the ratio of the output voltages at the detectors, and it is given by 








[%2)x(/l,x)R'2(/l)g(^)^.A 



(4.9) 



where e(A) is the emissivity, t(2, x) is atmospheric transmittance, and R'(2) is the photon 
detector responsivity. The subscript on the photon detector responsivity denotes the 
spectral band used in the measurement. In tliis equation, assuming that the wavelength 
parameters are spectrally constant and two identical detectors are used, the ratio be- 
comes identical to the ratio of photon emittance for black body (3.6). In order to gen- 
erate the output ratio calculation under more general conditions, the numerical method 
for integration known as Simpson's rule [Ref 12: pp. 95] is used. The algorithm for this 
calculation is shown in Appendix D. 
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1. Discussion of the computer simulations 

As described in this chapter, the actual ratio for the photon detector is the ratio 
of the output voltages. In the rest of this section, the simulation for the ratio calculation 
is discussed. Although the photon detector responsivity has wavelength dependence, as 
discussed in Section C, it is not a very sensitive factor up to the detector cut-off. 
Therefore, this is ignored in this discussion, and emissivity and atmospheric 
transmittance are considered. Since emissivity depends on the material of the thermal 
source, many kinds of emissivity values are possible. For this simulation, an arbitrary 
created emissivity curve shown in Figure 11 was used. For atmospheric transmittance, 
the atmospheric conditions of 1962 U.S. standard atmosphere were selected. The data, 
shown in Figure 10, are obtained from NFS LOWTRAN6 program. For the purpose of 
comparison with the results for ideal case, the same band pairs are selected as ones in 
Chapter III. The result for overlapping case is shown in Figure 12, for partially over- 
lapping case in Figure 13, and totally contained case in Figure 14. The curves generated 
in these figures correspond to 1 km and 5 km. With the exception of Figure 14, there 
is a one-to-one correspondence between measured ratio and temperature for each range. 
In Figure 14, there is a slight positive shape after T= 1300 K which violates the 
uniqueness condition. Although the values are shifted slightly, the shapes of the curves 
are essentially the same as ones for ideal ratio case shown in Chapter III. The general 
features of each case are the same. However, the results demonstrate that the curves 
depend on distance. This is quite inconvenient. The application of two-band temperature 
determination is then impossible without range information. Also, to lesser extent, the 
application is limited since the emissivity is material dependent. If the source for thermal 
radiation is not a gray body, then the emissivity of the source needs to be known in ad- 
vance. 
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Figure 10. Atmospheric transmittance for 1962 U.S. standard atmosphere: The 
path lengths are taken to be 1 km and 5 km in the horizontal direction. 
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Figure 11. Example of the emissivity (arbitrary chosen): This data does not cor- 
respond to the specific material. 
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Figure 12. 



The ratio of output voltages for nonoverlappiug case: Bands are (5.0 

nm, 4.0 /xm) and (4.0 fim, 3.0 nm). The distances are 1 km and 5 km. 
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Figure 13. The ratio of output voltages for partially overlapping case: Bands are 
(5.0 ^m, 3.5 ^m) and (4.5 ^m, 3.0 ^m). The distances are I km and 5 
km. 
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Figure 14. The ratio of output voltages for totally contained band case: Bands are 
(5.0 /rm, 3.0 ^m) and (4.5 /im, 3.5 /im). The distances are 1 km and 5 
km. 
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V. THERMAL RANGING TECHNIQUE 



From the previous chapter, the ratio of output voltages from dual spectral bands is 
not the same as the ratio of black body photon emittances under realistic and nonideal 
conditions. To solve and apply the method under these conditions can be quite difficult. 
On the other hand, one significant idea is hidden in this problem. 

The concept presented in this chapter is created from the difficulty, in the temper- 
ature determination method, produced by nonconstant atmospheric transmittance. 
Among assumptions discussed in Chapter IV, atmospheric transmittance breaks the va- 
lidity of temperature determination procedure because it depends on wavelength, the 
distance between an object and the detector, and atmospheric conditions. Although 
emissisivity and detector responsivity have to be taken into account, under some rea- 
sonable assumptions they can be neglected. It will be shown that the dependency of at- 
mospheric attenuation on the distance makes the concept of thermal ranging feasible. 

The use of infrared thermal signals has the advantage that active interrogation is not 
required since the object radiates energy or photons. This passive property of infrared 
detection can be significantly important in military situations. However, the passive de- 
tection of the target does not automatically provide range data and this is a disadvantage 
of thermal infrared detection systems. One interactive (nonautomatic) passive system 
involves the use of a trained human observer. The training allows the observer to rec- 
ognize a target and by noting the size of the target on the screen make an experienced 
estimate of the range. Currently ranging is done by active sensor systems such as radar 
or laser range finder. These active systems are capable of being sensed by the target. 
Although active ranging systems tend to provides precise distance information, they ex- 
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pose the instrument generating the signal. Therefore, the development of a passive 
ranging system would provide significant military advantage. 

The Electronic Support Measure (ESM) equipment in electronic warfare systems is 
the passive sensor in the microwave region. This gives target identification, direction of 
the target and also a very rough range information. For this method, the power trans- 
mitted by a radar of a target, the antenna gain, and attenuation factor must be obtained 
in advance. It seems that the same concept can be applied to infrared systems only if 
the temperature of a target, surface area, emissivity, and atmospheric transmittance are 
known. Therefore, this method is not practical. 

The follow'ing discussion proposes a new concept of passive ranging by infrared 
systems. It will be shown that this concept is a feasible solution to the problem of target 
ranging using only automatic passive methods. The following discussion consists of 
three steps: 

1. Step 1: Atmospheric attenuation coefficient is assumed to be constant. Temper- 
ature of the target is known. 

2. Step 2: Atmospheric attenuation coefficient is assumed to be constant. Temper- 
ature of the target is unknown. 

3. Step 3: Atmospheric attenuation coefficient is not constant. Temperature of the 
target is unknown. 

The concept of ranging shown here is based on the distance dependency of atmospheric 
transmittance. Step 1 and step 2 are intermediate procedures and they do not allow 
changing atmospheric conditions due to, for example, an elevated transmission path. 
In step 3, these prohibitions are rescinded and the ranging model has greater applica- 
bility. 

A. RANGING FOR KNOWN TEMPERATURE CONDITION WITH 
HOMOGENEOUS ATMOSPHERIC CONDITIONS 

In the analysis of this section the emissivity of the target and photon detector 
responsivity are assumed to be spectrally constant. In addition, it is assumed that at- 
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mospheric conditions are not changing along a line-of-sight trajectory between target 
and passive sensor system. Therefore, atmospheric attenuation coefficients will not de- 
pend on range but still exhibit wavelength dependence. By appropriate band selection, 
i.e., from a flat portion of the transmission curve [LOWTRAN], the wavelength de- 
pendence of the transmission can be approximately neglected with in each of the selected 
bands. Therefore, atmospheric transmittance for the /-th band is given by 

t/x) = exp( - trpc) (5.1) 

where a, is the atmospheric attenuation coefficient within the /-th band, and x is the 
transmission path length. Applying the assumptions stated to the ratio (4.9) leads to 

= (5.2^) 

After substitution of (5.1), the ratio can be written as 

= exp ( — (oj — 02 )x]^^ {5.2b) 

where is the ratio of the measured output voltage at the detector, is the ideal ratio 
of photon emittance, and a, are the atmospheric attenuation coefficient for the measured 
bands. In (5.2b), the ratio of the output voltages is obtained by direct measurement. The 
atmospheric attenuation coefficient is calculated from the atmospheric conditions, and 
the ratio of photon emittance is determined by the temperature given. Therefore, the 
transmission path length is the only unknown parameter. From (5.2b), the transmission 
path distance is obtained as 



X = 



1 

( 7 1 — <72 



In 






M 






(5.3) 
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The use of (5.3) makes passive thermal ranging feasible. However, it is applicable to very 
hmited conditions. In the battle field, the targets are usually unknown, and consequently 
the temperature of the target is hard to obtain. This means that the system should be 
capable of determining both range and temperature to be practical. The new proposed 
concept is described in the following sections. 

B. RANGING FOR UNKNOWN TEMPERATURE WITH HOMOGENEOUS 
ATMOSPHERIC CONDITIONS 

In this section a concept to determine the range and the temperature of the target 
is discussed. As in Section A, the atmospheric conditions are assumed to be homogene- 
ous in space. Mathematically, two nondegenerate equations can be used to solve for two 
unknowns. In this discussion, the two unknowns are range and temperature of a target. 
To get two equations, two distinct sets of ratio calculations are required. For that pur- 
pose, no less than three sets of band measurements ( i.e., band 1, band 2, and band 3 ) 
are required to produce two distinct equations such as (5.3). All the assumptions de- 
scribed in the previous section are also applied in this section except that the temper- 
ature of a target is unknown. The mathematical derivation to get the range equations 
is exactly the same as in the previous section. The two range equations for three sets of 
band measurements are 



X = 



1 

C7l — 02 



In 



0t\2 



(5.4) 



X 



1 

02 — 02 



In 



^23 



(5.5) 



where the subscripts serve as an index to identify spectral bands. Combining (5.4) and 
(5.5), one equation to solve the temperature is obtained as 
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(5.6) 
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a, -02 



In 






1 

02 -<73 
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However, this equation is impossible to solve analytically. A numerical method can be 
used to solve this equation. After getting the temperature, the range is easily obtained 
by using (5.4) or (5.5). Since the assumption that the atmospheric conditions are ho- 
mogeneous in space is only approximately valid for very short ranges, this method is not 
practical. The following section solves the difficulty of the determination of the range 
without assuming homogeneous atmospheric conditions. 



C. RANGING FOR UNKNOWN TEMPERATURE WITH INHOMOGENEOUS 
ATMOSPHERIC CONDITIONS 

For the following discussion, several equations are recalled for convenience. For 
four bands, the ideal ratios calculated from photon emittance are defined as: 



J _ Qi^\t ■^ 2 ’ ^ 
2 (/'- 3 .> 4 . 7 ) 



(5.7«) 



and 



,/ _ ■^6) T) 

Qih>h^T) ’ 



where the total photon emittance is given by 



(5.7i.) 



= T)dX, 



(5.8) 



and spectral photon emittance Qt can be expressed in terms of the gray body distribution 
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Inec 



(5.9) 






/[ cxpihclkT) - 1] ' 



With two band measurements, the temperature can be determined. With three band 
measurements, it should be possible to determine both temperature and range. As pre- 
viously discussed, range information is available as long as the effective atmospheric at- 
tenuations in each band are distinct. The problem posed by obtaining range information 
for a thermal source at the same elevation as the receiver (homogeneous atmospheric 
conditions) is not as complex as that associated with arbitrary elevation. Elevation 
changes in the trajectory of the radiation usually require consideration due to the cur- 
vature of the earth. In the general problem to be solved, both the receiver and the ther- 
mal source take arbitrary locations in the plane defined by the vertical line passing 
through the receiver and thermal source. This is represented in Figure 15. 

In the algorithm to be discussed, three arbitrary but distinct radiation bands are se- 
lected. Initially disregarding the effects due to the atmospheric attenuation, calculated 
ideal ratios of signals are given by 



Ideal signal Band 1 
= Ideal signal Band 2 ’ 



(5.10a) 



I Ideal signal Band 3 

^3i(^ Ideal signal Band 1 



( 5 . 10 ^ 7 ) 



Figure 16 shows the temperature dependence of these curves for bands defined by: 

(.Ij, > 12 ) = (12.0 fim, 10.8 fim) 



(y? 3 , /I 4 ) = (10.8 nm, 9.6 uni) 
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(-^s. = (9-6 M"*' 8-2 fim) 



The atmospheric attenuation in the form of transmittance is provided from the NPS 
LOWTRAN6 program (Ref 13] using data corresponding to a vertical path. Atten- 
uation coefTicients are derived from 



where i= 1,2,3 indicate the wavelength band and j is the index of the specific layer. The 
quantity is the length of the trajectory in the y-th atmospheric layer. See Figure 17. 
Each fj has the same length if the line-of-sight angle is 90 degree. The attenuation coef- 
ficient o,y is obtained by taking the 90 degree angle. The standard band average for Xy 
is given by 



With reference to Figure 18, the effective approximate transmission factor t„ for the i-th 
band, is given by: 



where is the number of layers between source and receiver, and m is the layer num- 
ber of the receiver. Here, r* is the length of the trajectory in the layer k + m — The 
approximation involved in (5.12) will be discussed at the end of this section. Expression 
(5.13) can be recast into the form: 



= Inepers/m] 



(5.11) 




(5.12) 




(5.13) 
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T; = exp( — oiR) 



(5.14) 



where 



Oi = 



^«x 

I 

A=1 






R 



(5.15) 



and when the range R must satisfy; 






ax 




k=\ 



(5.16) 



With the exception of the first and final layer, the value for each r, is obtained from the 
rule 



n = 




(5.17) 



as showm in Figure 19. Each R, is calculated from the Cosine Law relation derived from 
the geometry of Figure 20 

(//, + R^f = Rf + (h, + Rsf - 2Rih, + R^) cosid + f ) (5.18a) 

or equivalently 

(Hi + Re)^ = r} + {h^ + R^f + 2Rlh^ + Re) sin 6, (5.18^) 
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where Q is the line of sight angle between receiver and source measured with respect to 
the horizontal. R, can be obtained from (5.18) using the quadratic formula with the 
correctly chosen ' + ' sign, f?, is also obtained from (5.18) after taking 

( 5 . 19 ) 

where the int[J results in the integer part. This follows from examination of the geom- 
etry in Figure 15. For /? 2 ,. /?* values for H* follow 

//ft = //i -t- (A - l)Ar (5.20) 

which are used in (5.18) to calculate subsequent /?*. For clarity, it is noted that: 

/?! = rj (5.21a) 

/?ft = Tj -f -f-.... + Tft (5.21Z>) 

until k = as defined by (5.16). Equations (5.17) combined with (5.20) are solved it- 
eratively to generate the values for r, except for Lastly, rf,„, is calculated from: 

Once r, are determined for all i up to , the effective transmission for each band (5.14) 
can be calculated from (5.15) and (5.16). The measured ratios, i.e., including the effects 
of the atmosphere, are related to ideal ratios (5.7) according to: 

= ( 5 - 23 «) 



and 
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(5.236) 



— ^31 



h_ 

’ 



in agreement with (5.2a). After substitution from (5.14), (5.23) can be written as 

= ^12 exp [ - (oi - <72)i?] (5.24a) 



and 



< = ^;,exp[-(a3-<r0i?]. 



(5.246) 



w'here R is the distance between thermal source and receiver. Figure 21 shows the 

measured ratios calculated from (5.24) for conditions: 

• <7, = 0.0001494 [neperslm] 

• <72 = 0.0001027 [neperslm] 

• = 0.0004622 [neperslm] 

• R=100 [m] 

Solving for R from (5.24), it follows that 



12 1 ^12 
r^\T) = In — , 



(5.25a) 



and 
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mxiT) 



(5.256) 



W’here the subscripts are employed to denote spectral bands. In the algorithm, it will be 
assumed that the correct measured ratios, consistent with (5.24), are provided by exper- 
imental testing. Ideal ratio curves can always be calculated according to the methods 
(5.10) previously discussed. Constants <7„ aj, <7j can be approximately calculated from the 
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NPS L0WTRAN6 data. By checking for the crossing point between the two curves 
(5.25), both the temperature, T, and the range, R, are uniquely determined as shown in 
Figure 22. 

The algorithm sketched out is not yet complete. The values of <t„ and ctj are de- 
pendent on the range, R, since the weights for evaluating the coefficients a, (5.15) depend 
on all the values for r,. In order to calculate the coefficients a, (5.15), an initial guess is 
made for a range candidate, within limits dictated by the domain of LOWTRAN 
data sampled. This value for /?co ts used in the calculation of <r„ and Oy The search 
for a numerical crossing point for curves (5.25) produces a new value for R which can 
be denoted as R^y The value for Rg, generated by this first iteration is expected to be 
closer to the actual range but with an unacceptable level of error. The process is re- 
peated by iterating with the R^y The next estimate generated by the algorithm, is 
tested against the previous one, R^y In general, the procedure is stopped when 



^Cm ~ ^ 



Gm—\ 



R 



Gin 



<5. 



(5.26) 



where 6 is the predetermined convergence parameter and the index m specifies the 
number of iteration cycles. In Appendix E, the algorithm for the process is presented. 

For the purpose of checking the algorithm, the range candidate, is selected as 
shown in Figure 23. Actual ranges are chosen from 10 meters to 700 meters and Rco’s 
are 700 meters for up to 350 meters of actual ranges and 10 meters for longer than 350 
meters. This choice appears to correspond to worst-case guessing. The results for error 
calculation are obtained for 6= 0 degrees, 30 degrees, and 45 degrees at a detector 
height of 5 meters. The actual temperature and convergence parameter b are chosen to 
be 353.54 K and 0.001, respectively. The error of range determination is shown in Figure 
24, and the error of temperature determination is shown in Figure 25. Both errors in- 
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crease as the actual range increases. However, the maximum error in temperature and 
range is insignificantly small. The maximum error in ranging is 0.09 % and the maximum 
error in temperature is 0.0008 %. This results shows that the method performs remark- 
ably well for the conditions tested. 

In the calculation of the transmission factor for the /-th band, the precise t,* should 
be calculated according to the weighted average 



rather than the standard average (5.12). It should be noted that, the flatter the photon 
emittance distribution within the spectral bands, the smaller the discrepancy between 
calculations (5.12) and (5.27). 

Despite the simplifying assumptions used in the model, the results demonstrate the 
feasibility of extracting both target range and temperature from thermal radiation 
measured over three bands. 




(5.27) 
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Figure 15. Deflnition of vertical line from the earth through the receiver 
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Figure 16 . Ideal ratio curves 
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Figure 17. Definition of layers in atmosphere 
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Figure 20. 




Geometry of ranging (cosine law) 
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Figure 21. Measured ratio curves 
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Figure 22. 



Graphical solution of range and temperature 
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Figure 23. Initial guess range (the worst case) 
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Figure 24. Error of range determination 
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Figure 25. Error of temperature determination 
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VI. CONCLUSIONS 



Previous work on dual spectral band measurements of radiant emittance demon- 
strated that, under somewhat ideal conditions, the ratio of the measurements could be 
used to determine the temperature. The main stipulated condition was that the bands 
were not self-contained. In the Chapter III of this thesis, the same conclusions were 
obtained for measurements of photon emittance. 

A general model encompassing arbitrary wavelength-dependent emissivity, 
responsivity, and transmission factor was proposed in Chapter IV. Although the general 
problem is not analytically tractable, the methods applicable to constant emissivity, 
responsivity, and transmission factor (Chapter III) are approximately valid if the spec- 
tral bands are narrow enough. However, from the practical view point of making a 
signal measurement in the presence of noise, the spectral bands should not be too nar- 
row. Therefore, the more general model is often needed under practical conditions. For 
reason explained in Chapter IV, the photon detector responsivity introduced is fairly 
insensitive to wavelength up to the cut-oIT wavelength. In addition, for a broad class 
of materials emissivity curves are quite flat. Therefore, the focus of the reminder of the 
thesis concentrated on effects due to atmospheric attenuation. Specifically, knowledge 
of the range to the target would be necessary to apply the method discussed in Chapter 
IV. 

The concept presented and tested in Chapter V is created from the requirement that 
the dual spectral band method of temperature determination would require range infor- 
mation to determine the temperature. This assumes a disparity in attenuation effects 
between bands. By employing three bands instead of two, this problem is circumvented. 
Additional complications created by changing atmospheric conditions with trajectory 
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path are handled within the model. The numerical atmospheric attenuation data used in 
the model testing was taken from LOWTRAN. The results for this preliminary model 
are extremely encouraging. Future work needs to be performed to develop and test a 
more rigorous model. 
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APPENDIX A. DERIVATIVE OF 01 WITH RESPECT TO T 



The purpose of this appendix is to prove equation (3.7) which appears in the main 
text. Equation (3.7) had been introduced in [Ref. 8: pp. 1257] without proof. Since the 
ratio of photon emittance is a function of a, it requires a few steps to get its derivative 
with respect T. For the reader's convenience, relations (3.2a), (3.5) and (3.6) are re- 
produced, which show that 






CaX2 






CiXXt 

J 






(A\) 



where 






exp{\p) - 1 



(A2) 



and 



a 







Following the product rule, the derivative of the ratio of photon emittance (Al) with 
respect to temperature is 
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8^ 

8T 



a f<«2 f“^2 A 

-^J gii'W-j si'Pyii/ 

a 0UC3 a 0UC3 



«3 



(^4) 



where the temperature dependence is expressed through a. One approach to simplifying 
the expression is to transfer the dependence of the integral limit on a evident in (A4) to 
the integrand. A substitution, 



ij/ = ax 



(ASa) 



will transfer this dependence. It follows that 



[ ^g(ax)dx. (ASb) 

aXi Xj 



Therefore according to the product rule, 



a Cox, fXj 

g{'l'W=j ^(axy.v + a -^g(ax)dx. 

OLX{ Xj Jy 



(A6) 



Now from (A2) and the transformation rule (A5a), 

^^"^)=expM-T- 



iAl) 



and it follows from (A7) that 
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a , . _ 2ax^ exp(our) 

6a ^ exp(ax) - 1 [ exp(ag-) - ’ 



Direct substitution of (A8) into (A6) leads to : 



{AS) 



J 'g{'P)d^ = f 'g{ax)dx ■ 

•'ax, 



rx, 



2{oixf 
exp(gjf) — 1 



C^J 



dx — 



(gx)^ exp(gx) 
C exp(gx)- 



■dx. (^9) 



Noting the equivalence between ^(gx) defined in (A7), and the integrand of the second 
integral (A9), it follows that : 



_6_ 

da 



J ^g{^)d^ = sj ^g{ax)dx - 

CLXt X: 



rx: 



(gx)^ exp(gx) 
C exp(gx) - 1 



■dx. 



{A\0) 



After reversing the transformation (A5a), the integrands in (A 10) will depend on \J/ and 
(AIO) leads to: 



da 



Caxj 3 Caxj , 

g{^l>)dyl,=^\ g{m—^ 



'P exp(^) 

[ exp(^)- 1]^ 



d\j/. 



(^ 11 ) 



For compactness of representation, define 



A^) = 



exp(<A) 

C expiij/)- 1]^ 



{A\2) 
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Then combining the result (A3) with (All) leads to : 



OX/ 



y[J ~ 

|_*'ooc, -^ax, 



iAl3) 



After substituting (A 1 3) into (A4) 



e^p 

dT 



[j y(i/')#-3f ^^(^)</ia 1| g(^)#- f V(^)# f yi^A)# - 3 f*^g(^^)#l 

1 L*^a J^olX 3 a J 



[i~ 



g{lp)d^ 



[A14) 



which after simplification leads directly to (A 16) or (3.7), as given in the text. 



d^r 



raX2 rax^ f«2 roa;, 

A^)d'lf\ gWd^-\ g^dipl A^)d'P 

1 *’a *’oU3 *’a **0X3 



dT T 


r f“^ 






gWd^ 


2 




L*’«3 J 





(^15) 
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APPENDIX D. ASYMPTOTIC APPROXIMATION FOR RATIO 



CALCULATION 



The purpose of this appendix is to prove equations (3.22) and (3.23) which appear 
in the main text where the asymptotic approximations are applied. 

A. FOR LARGE ALPHA 

For large a (i.e., small temperature), since x, values are fixed, \j/ takes on large posi- 
tive values Therefore, and functions are approximated from (3.5) and (3.8) as 
follows: 




(51) 



and 




(52) 



After direct substitution into (3.8), it follows that; 




a 



3 2 

4/ exp( — exp( — 

0X3 




One round of integration by parts yields 
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/'oUCj 



exp(-^)]r^ + 3 



exp( - ^)di> 



*0X4 



\J/^ exp( — ^)d\}/ 



C0LX7 



exp( - exp( - + 3 



«3 

/• 0 OC 4 



• m 



V exp( - 



<«3 



After algebraic simplification, (B4) becomes 



lf,g]^[. - exp( - 



•axt 



exp( — ^)d\]/ 



«3 



f^OLX'i 



\}/^ exp( - i/)d\p\_ - \}/^ exp( - 



(55) 



According to [Ref. 12: pp. 85] 



exp( — \(/)d\l/ = — \{P' exp( — ^) — 2\p exp( — \j/) — 2 exp( — \p), (56) 



which, for large positive ^ (B6), is approximately given by 



2 2 
\l/ exp( — \l/)d\j/ii — ^ exp( — iff). 



{Bl) 



After substitution of the approximate form (B7) into the integrals (B5), [f,g] becomes 
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lf,g]-U'^ exp( - exp( - - [<A^ exp( - exp( - (A)]^-(^ 8 ) 



which, after simplification, becomes 




Since, for large a, the second terms within each bracket of (B9) become dominant, [f,g] 
is approximated as 



This last form is applied in Section C of Chapter III. 

B. FOR SMALL ALPHA 

For small a (i.e., large temperature ), ^ is small. Thus, g(^) and/(^) are approxi- 
mated from (3.5) and (3.8) by binomial expansion. After application of binomial expan- 
sion, g{ij/) is integrated as 



[/;g]=:a^[jC 3 exp( - ax^) exp( - a) -x^ exp( - ajC 3 ) exp( - a)]. (BIO) 



Recombination of (Bll) leads to: 



[/;g]=ia^X 3 (l - JC 3 ) exp [ - a(x 3 4-1)]. 



(Bll) 




(512) 



And a binomial approximation for small ip leads to: 




(513) 
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which is integrated as 






(514) 



Similarly, after applying binomial expansion to (3.8), the integration of J[\l/) becomes 






(515) 



Algebraic expansion of(B15) leads to: 



J/t'/'W 




1 

2 



3 ,2 , 1 ,3 

"I?'*' +lo 



)• 



(516) 



which is approximated as 






(517) 



Therefore, for small a, [f,g] becomes 




(518) 



which is expanded as 
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LAg ]- -^2 (-^2 - -^ a ) - ^ aVi - 1)(^4 - -^ a ) “ ^ (-^2 - ^ 3 ] 



32 



32 



- -^ (2^4 - ^a)(^2 - 1) + 
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6/4 4w 2 ,N ft / 4 4w 3 
a (^4 -•^ a )(^2 “*) “ U (^4 - ^ a )(^2 “ 



Since, for small alpha, oc* terms are dominant, [f,g] approximated as 

l/'.^l- -J2 - ^a)(^2 -1) - (^2 - 1)(^4 - 4)1 



This final form is applied in Section C of Chapter III. 



■ ( 519 ) 

1 ) 



.. (520) 
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APPENDIX C. ALGORITHM TO CALCULATE THE RATIO OF 

PHOTON EMITTANCE 




photon or radant 
omittance dstribution 
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APPENDIX D. ALGORITHM TO CALCULATE THE RATIO OF 
OUTPUT VOLTAGE AT THE DETECTOR 
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APPENDIX E. MODEL FOR PASSIVE RANGE AND TEMPERATURE 

DETERMINATION 



A MAIN ALGORITHM 
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B. LOWTRAN DATA GENERATION 



Inputs 
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C. LAYER (SUBROUTINE) 



Inputs 
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D. MEASURED RATIO PREDICTOR 
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